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^ Abstract 

We examine a possibility that, when a black hole is formed, the 
Q^ I information on the collapsed star is stored as the entanglement entropy 

' between the outside and the thin region (of the order of the Planck 

length) of the inside the horizon. For this reason, we call this as the 
Q^ I entanglement entropy of the black hole "horizon" . We construct two 

models, one is in the Minkowski spacetime and the other is in the 



O . Rindler wedge. To calculate the entropy explicitly, we assume that 

I , the thin regions of the order of the Planck length of the outside and 

^j-jl inside the horizon are completely entangled by quantum effects. We 

also use a property of the entanglement entropy that it is symmetric 
under an interchange of the observed and unobserved subsystems. Our 
r> I setting and this symmetric property substantially reduce the needed 

■ numerical calculation. As a result of our analysis, we can explain 

the Bekenstein-Hawking entropy itself (rather than its correction by 
matter fields) in the context of the entanglement entropy. 



1 Introduction 

There is a well-known analogy between black hole physics and thermody- 
namics. This fact is called the black hole thermodynamics [|l| . In particular, 
as first pointed out by Bekenstein 0, we can think of the area of the black 
hole horizon as the entropy (up to a proportional constant) by using the area 
theorem which states that the area of the black hole horizon does not 
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decrease. Since the black hole emits thermal radiation of matter, which is 
called the Hawking radiation |, ^ , we can decide the temperature of the 
black hole. Thus, the entropy of the black hole is calculated as 



^=I7^A (1-1) 

where A is the area of the horizon and Ipi = {hG/c^Y^"^ is the Planck length. 
This is the Bekenstein-Hawking entropy. 

There have been many attempts to understand the origin of this black 
hole entropy: For example, those considerations on the basis of the value 
of the Euclidean action 0, ||, ^ [l^, |ll|, the rate of the pair creation of 



black holes |12|, the Noether charge of the bifurcate Killing horizon [13 



IJ] or the central charge of the Virasoro algebra |]T5|, 0, |T3- Among past 



considerations, we consider the entanglement entropy |18, 19, EG, 21, 22, 23 



2J, |25|, |2^ as the most attractive candidate for the black hole entropy. The 
entanglement entropy is the measure of the information loss due to a division 
of the system; this direct connection of the entropy with the information loss 
is not clear in some other approaches to the black hole entropy. If we divide 
the system into two subsystems A and B, and ignore the information about 
B and observe only A, we can view the pure state of the total system as an 
effective mixed state for the subsystem A. The entanglement entropy is the 
von Neumann entropy of this effective mixed state. If the original pure state 
is an entangled state, the entanglement entropy is non-zero. On the other 
hand, if the original pure state is not an entangled state, the entanglement 
entropy is zero. That is, if the original pure state is not entangled, there is 
no information loss when we ignore B. Note that the entangled state and 
the entanglement entropy is a purely quantum mechanical notion and there 
is no counterpart in classical physics. 

When the concept of the entanglement entropy is applied to the black 
hole, it measures the information loss due to a spatial separation. Most 
previous works on the entanglement entropy were concentrated on the en- 
tanglement between the bulk regions outside and inside of the black hole 
horizon. In this paper, we instead discuss the entanglement between the 
outside and a thin region (of the order of the Planck length) inside of the 
horizon. For this reason, we call this as the entanglement entropy of the 
black hole "horizon" . 
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We consider that this approach is justified physically by the following 
discussion: For simplicity, we consider a quantum field on the extended 
Schwarzschild spacetime rather than a dynamical spacetime which describes 
the gravitational collapse to the black hole. Since we want to calculate the 
entropy of the black hole itself, the quantum state of the field must be a 
"vacuum". We thus consider the Killing vacuum, which is defined by using 
the Killing time. (Note that the Kruskal vacuum, which is defined by us- 
ing the Kruskal time, contains the thermal radiation of the Killing particles, 
namely, the Hawking radiation. Therefore, if we chose the Kruskal vacuum, 
the resultant entropy would be considered as the entropy of the black hole 
and its correction by the matter field.) Since the Killing vacuum is expressed 
as the tensor product of the states in one asymptotically fiat region I (out- 
side) and the other asymptotically fiat region II (inside), the entanglement 
entropy between the inside and outside of the horizon becomes zero. 

However, if we consider the effect of the quantum gravity, this vanishing 
entanglement entropy is not true any more since we can not divide the system 
sharply due to the quantum fiuctuation of the horizon. By this correction, 
the Killing vacuum is deformed to some entangled states between the inside 
and outside of the (classical) horizon. The depth of the entanglement is of 
the order of the Planck length because this is the effect of the quantum grav- 
ity. Of course, it is difficult to achieve such a calculation. To estimate this 
entanglement entropy, we first assume that the thin regions of the order of 
the Planck length outside and inside the horizon are completely entangled by 
quantum effects. Namely, the main features of the states in the thin region 
inside the horizon are smoothly extrapolated from those of the outside the 
horizon. The major ansatz of our calculation is that the entanglement en- 
tropy between the thin region inside the horizon and all the states outside the 
horizon is approximated by the entanglement entropy of the Killing vacuum 
between the thin region of the order of the Planck length outside the horizon 
and the rest of the states outside the horizon. That is, we consider that the 
effect of the quantum gravity is approximated by the shift of the (classical) 
division of the system rather than the deformation of the state. The shift 
is of the order of the Planck length because this is induced by the quantum 
fluctuation of the horizon. 

The present ansatz is analogous to the setting which has been considered 
in Ref.[^ in a different context. They have considered a thin spherical shell 
infalling toward a Schwarzschild black hole and the entanglement entropy 
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of the Killing vacuum associated with the division by a timelike surface 
which becomes the horizon after the passage of the shell, but it is in the 
Schwarzschild spacetime before the passage of the shell. Our case is differ- 
ent since we analyze the entanglement entropy generated by quantum effects 
after the formation of the classical horizon. Also, the calculation becomes 
much simpler in this paper. The key point is that the entanglement entropy 
is symmetric under an interchange of the role of the subsystems A and B. 
Our setting and this symmetric property make the calculation very simple 
and substantially reduce the needed numerical calculation. Moreover, since 
the calculation in this paper is based on the Bombelli-Koul-Lee-Sorkin type 
calculation |jl8| rather than the Srednicki type calculation |]19|, we can find 
directly that the entanglement entropy is proportional to the area without 
plotting the entanglement entropy to the area. 

There are some comments on the above setting of the calculation in this 
paper: As is well known, Euclidean geometry plays an important role in the 
Gibbons- Hawking method || or some other Euclidean approaches to the 
black hole entropy. Especially, the temperature of the black hole can be well 



understood as the period of the Euchdean time in Euclidean geometry |27, 28 



29| . Since the entropy is the conjugate variable to the temperature, we want 
to understand it within Euclidean geometry. However, the Euclidean black 
hole does not have the "inside" of the horizon p, 0. On the other hand, 
in Euclidean gravity, the horizon is the fixed point of the Euclidean time 
translation, called the bolt and an obstruction to the foliation by the 
Euclidean time. Therefore, to achieve the Hamiltonian formulation, we want 
to eliminate the degrees of freedom near the horizon |3^, For this reason, 
the above setting of our calculation appears to be reasonable if we persist 
on the Euclidean picture. Then, we notice that the energy, temperature and 
entropy could be understood in relation to the Euclidean time translation: 
That is, we regard that the energy is its charge, the temperature is its period 
and the entropy is concerned with its fixed point. 

Moreover, to reach the horizon, we need an infinite time if the "time" 
is measured by the asymptotically Minkowski time (not the proper time). 
Thus, we can consider that the horizon is a "boundary", at least, for the 
observer at infinity. However, we do not impose any boundary condition at 
the horizon since we do not make any measurement there. Consequently, it 
is natural to take the summation over the state at the horizon [Q. 

The plan of this paper is as follows. In Sec.^ we briefly review the notion 
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and basic properties of the entanglement entropy, and then derive a basic 
formula to calculate it. In Sec.|, we construct two models and calculate the 
entanglement entropy explicitly. In Sec.^ we conclude and discuss the results 
of this paper. 

2 Entanglement Entropy 

We review the notion and properties of the entanglement entropy and then 
derive a basic formula |jT8[ to calculate it. 

2.1 Definition 

Let us consider the case where the total system can be divided into two 
subsystems. Then, the Hilbert space of the total system Ti can be written 
by the tensor product, 

U = Hi®n2. (2.1) 
A state |\E') G 7i is called entangled if the state can not be written as 

l^) = 1^1)1^2), (2.2) 
where G Tii and \ip2) G 7^2- For example, if |a), \h) G Tii and |q;), \(3) G 

m = \a)\a) + \bm (2.3) 

is an entangled state and 

1^) = \a)\a) + 2\a)\l3) + \h)\a) + 2\h)\p) 

= (|a) + |6))(|a) + 2|/?)) (2.4) 

is not an entangled state. 

Moreover, we assume that we are going to ignore the degrees of freedom 
of 0,2- To achieve this, we define a reduced density matrix pred for Tii from 
the (pure) state of the total system |\E'), whose matrix elements are given by 

(a|p,ed|6) = E(H(«l)l^)(^l(l (2-5) 
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where \a),\b) are the arbitrary states of Tii and {|a)} are the orthonormal 
basis of 7^2- Then, the expectation value of an operator O which acts only 
on 7^1 becomes 

(^|C>|^) = Tri(predO), (2.6) 

where the trace is taken over the states of Tii. By this way, as far as the 
subsystem Hi is concerned, the pure state of the total system can be 
viewed as the mixed state pred- 

Now, the entanglement entropy is defined by the von Neumann entropy 
of this reduced density matrix, 

5'l2 = -Tri(pi.edlnpred) 

= -"^Pn^iapn, (2.7) 

n 

where {pn} are the eigenvalues of p^ed- Note that the range of the entangle- 
ment entropy is 

0<Si2<\nN, (2.8) 

where N is the dimension of Tii. 

If the original state |\E') is not entangled, pred remains pure and, thus, 5*12 
becomes zero. On the other hand, if |\E') is entangled, pred becomes a mixed 
state and 5*12 is nonzero. Thus, the entanglement entropy is a measure of the 
entangled nature (or EPR correlation) of the original state. 

For example, let us consider the system which consists of two spin-1/2 
particles: 

If a state of the system is an EPR state, 

l^) = ^(ITi)IT2) + Ui)U2)), (2.9) 

then the reduced density matrix becomes 

Pred=(^ l), (2.10) 

and the entanglement entropy is S12 = In 2. This state has the maximum 
entanglement entropy |33|. Note that, since there is a perfect EPR correlation 



between these particles, we can get full information about one particle by 
an observation of the other particle. Thus, these particles are maximally 
entangled. 
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On the other hand, if a state of the system is not an entangled state, 



^(ITi) + lii))®^(IT2) + U2)) 

^|Tl)|T2) + |Tl)|i2) + |il)|T2) + |il)|i2)), (2.11) 



2 

then the reduced density matrix becomes 



Pred = ? ? , (2.12) 
^2 2' 

and the entanglement entropy is S\2 = 0. This state does not have the 
entanglement entropy. Note that, since there is no EPR correlation between 
these particles, we can not get any information about one particle by an 
observation of the other particle. Thus, these particles are not entangled. 

Note that one of the important properties of the entanglement entropy is 
that it is symmetric under an interchange of the role of Tii and 7^2, 

5l2 = ^21. (2.13) 

This is because the entanglement entropy measures the EPR "correlation" 
between two subsystems, which is symmetric by definition. As for more 
detailed analysis, see Refs. |jl9|, ^ . 



2.2 Basic Formula 

Let us consider a system which consists of coupled oscillators, {q^^- Now, we 
will calculate the entanglement entropy of the ground state when the system 
is divided into two subsystems, {(f^ and {g"} flSf . 



The Lagrangian of the total system is given by 

L = ^Gab q^q"" - \vab g^g^. (2.14) 

(We assume that Gab and Vab are symmetric and positive definite matrices 
of constants.) The canonical momentum conjugate to q^ is 

PA = GABq''. (2.15) 
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By using (G ^) which is the inverse matrix of Gab defined by 



{G-y^Gsc = S%, 



the Hamiltonian becomes 



H^\iG-Y^PAPB + \vABq^q^. 
Moreover, we define Wab by 

(G-Y^WacWbd = VcD. 



(2.16) 



(2.17) 



(2.18) 



That is, Wab is the square root of Vab in terms of the metric {G ^)^^. Then, 
by using the canonical commutation relation 



Bi 



(2.19) 



one finds that the Hamiltonian becomes 
H = \{G-Y'' {pa + iWAcq'') {PB - iWEDq"") + ^iG-yWAB- (2.20) 

Thus, we can define the creation operator and the annihilation operator 

ua by 



OA = {pa - iWAcq^) , 

= -^{pA + iWAcq^). 
The commutation relation between these operators are 



We can then write the Hamiltonian as 



^Wab. 



1\AB 



+ 2^AB 



(2.21) 

(2.22) 

(2.23) 
(2.24) 



The first term is the number operator and the second term is the zero-point 
energy. 
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The ground state of this system is given by 

1 



V2 



PA-iWAcq"-] |0) = 



(2.25) 



The wave function of the ground state is obtained by 

/ d 



+ w^AcgM (K}|o) = o, 



(2.26) 



since pa = —id/dq^ in the Schrodinger representation. The normahzed 
solution of this equation is 



({/}|0) 



det 



W 



TC 



1/4 



exp 



The density matrix of this ground state is 
piiQ^'yAQ"'}) = ({/}|0)(0|{O) 



det 



W 



1/2 



exp 



(2.27) 

(2.28) 
. (2.29) 



/^^9>({9^9"},{9'^0) (2.30) 



Now, we divide the system {q^} into two subsystems, {g"} and {q"'}- If 
we want to ignore the information on {g"}, we take the trace over {g"} and 
consider the reduced density matrix as Eg. (|2.5|), 



Prod {{q''},W'} 
By dividing Wab into four blocks 

Wab= ( 

we find that 



(2.31) 



det — 

vr . 



1/2 



exp 



X exp 



(2.32) 
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where 



Mab = 
Nab = 

and we have used that 

detiy = det 

= det 



Aab 

(bd-^b 



BD-^B'^ 



ab 



(2.33) 
(2.34) 



ab 



{B^)^b D 



A, 



ab 



BD-^B^ 




ab 



B, 

D. 



afi 



detMdetD. 



Moreover, we can choose a basis {q""} in which both Mab and Nab are 
diagonaL Then, in this basis, the reduced density matrix becomes 



Pre, {m, w'}) = nj;^ [-1 ((^~")' + 

where {Xa} are the eigenvalues of the operator 

A\ = {M-^Y'N,b. 



(2.35) 



(2.36) 



In order to obtain a simpler expression for A'^^, we divide the inverse 
matrix of total Wab into four blocks. 



-l\AB 



j^ab j^aP 
^^Tyb jjap 



By definition. 



A"'^Bb^ + B^f^Dp^ = 
From Eg. (|2.39|) , we can find 

j^ap ^ -A^^Bb^iD-y^. 



(2.37) 



(2.38) 
(2.39) 



(2.40) 
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Combining this with Eg. ( |2.38| ), we obtain that 

(M-i)"^ = (2.41) 

Then, from Eg. ( |2.40| ) and Eg. ( p.41| ), it is easy to see that Eg.( |2.36| ) becomes 

A'^, = = A'^'A,, - 5\. (2.42) 

Note that the total reduced density matrix can be written by the tensor 
product, 

Pred = (g)P0(Aa), (2.43) 



where 



Po(A) 



exp 



Thus, the entropy is given by the summation with respect to each Aa, 



S = -Trpred Inpred = 'S'(Aa) 



where 



5(A) = -Trpo(A)lnpo(A). 
To calculate S{X), we must obtain the eigenvalues of po(A), 



dq' Po(A; q, q')fn{q') = Pnfn{q)- 
By using the formula for Hermite polynomials p4 |. 



r dx e-("-^)'if„(ax) = (l - a^Y^^ 

J —oo ^ ' 



ay 



.(1 -a')'^' 

we find that the eigenvalues and eigenfunctions are given by W^ 



Pn 

fn{q) 



[1 - p)p^ 



exp 



-^iq 



HniViq), 



where 
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Vl + A, 



A 



(2.44) 

(2.45) 
(2.46) 
(2.47) 

(2.48) 

(2.49) 
(2.50) 

(2.51) 



Then, the entropy for A is given by 



5(A) = - lnp„ = - ln(l - /i) - In/i. (2.53) 

In summary, in order to calculate the entanglement entropy for the ground 
state of coupled oscillators, one must first obtain the eigenvalues {Aa} of 
in Eq.( |2.42| ). Next, for each eigenvalue A^, one has to calculate Ha by 
Eq.( |2.52|) . Finally, the entanglement entropy is given by 



5 = E 



- ln(l - /ia) - , ^° In Ha 

1 - Ata 



(2.54) 



3 Models 



In this section, we will construct specific models and calculate the entan- 
glement entropy. We consider a free scalar field in a background spacetime. 
Since the field can be viewed as a set of coupled oscillators, we can use the 
formula in the previous section. 

We must divide the set of oscillators into two subsets. In most previous 
works, it was divided into the oscillators outside and inside of the black 
hole. Instead, in this paper, we will divide the system into the oscillators 
completely outside and within a thin region A around horizon, based on 
the discussion in Sec.|l[ Since the thin region is induced by the quantum 
fluctuation of the horizon, the width of the region A is of the order of the 
Planck length. Namely, a ^ Ipi, where a is the width of the region. 

Furthermore, we can make the calculation simpler. If one applies the 
conventional calculational scheme to our setting, we ignore the degrees of 
freedom inside of A (near the horizon). Instead, we here ignore the degrees 
of freedom outside of A in this paper. Of course, this gives the same entangle- 
ment entropy in our setting, because the entanglement entropy is symmetric 



as in Eq. (|2.13| ). Moreover, since the width of the region is of the order of 
the Planck length, we can treat the field within A as a single oscillator if 
we assume a momentum cut-off associated with a, which is of the order of 
the Planck scale. (If we assume a different momentum cut-off, which is still 
of the order of the Planck scale, we must consider the field within the thin 
region as a set of oscillators rather than a single oscillator. However, the 
number of the oscillators are still of the order of 1. Although this would 
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change the numerical value of the coefficient of the entanglement entropy, it 
would be still the same order. Thus, the final conclusion is unaffected.) This 
makes the calculation quite simple. Especially, the matrix A"^ becomes 1x1 
matrix and the eigenvalue Aq, itself. 



3.1 Simple Model 



First, we consider a free scalar field in the flat spacetime and adopt the 
Minkowski coordinates {t,x,y,z) = {t,x). The index A is now replaced by 
X. We assume that the "horizon" is at x = —L (— — oo). Then, the thin 
region near the horizon becomes A = {{x,y, z)\ — L < x < —L + a}. (See 
Fig.|I|.) Since the action is 



S 



dtd^x — 



{dt(t)f - {d(l)f - m^(P^ 



(3.1) 



one can easily find that Gab and Vab appearing in Eq.( p^.l4D become 



G{x, x') 
V{x, x') 

Then, one finds that 

W{x,x') 

W-^{x,x') 

Thus, k\ in Eq.(|2^) becomes 



G ^{x, x) = S{x — x') 



(2vr)^ 



{\k\' + 



J (27r)3 V ' 

kr + 



"'■/^ ik-{x—x') 



A(x, x') 



dx" 



W-\x,x")W{x",x') -5{x-x'), 



where x,x' G A. 

To solve the eigenvalue equation. 



/ dx'A{x,x')F{x') = XF{x), 

J A 



(3.2) 
(3.3) 



(3.4) 
(3.5) 

(3.6) 



(3.7) 
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we make the ansatz 

Fix) = e^P-^'fix), (3.8) 
where x = {y, z) and p = {py,Pz)- Then, the eigenvalue equation reduces to 

-L J-L J 277 J 2n 

=^—^K' + M2 e^'=^(--") e^K(^"-^') f(x') = (A + l)/(a;), (3.9) 
where 

Mp = ^m? + |p|2. (3.10) 
Moreover, we use an approximation, 

J 1^^" dxG{x) r-^ aG{-L + a/2), (3.11) 

which corresponds to the prescription that we treat the field within A as a 
single oscillator. We then find that 



\ + l = a' , Jk'^ + M2. (3.12) 



This would diverge unless we introduce a momentum cut-off kc- The mo- 
mentum cut-off can be decided by 

J-L J-oo 27r 



~ a 



dkx 

-kc 271 



(3.13) 



5 

71 



in relation to the width of the region and the approximation Eq.( |3.11| ). Thus, 

kc = -. (3.14) 
a 
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By using this momentum cut-off, we obtain that 



A(C) 



In 



v^r+c-i 



-1, 



where ( = Mpa/n. Then, from Eq.( p.52|) and Eq.( p.53|) , one finds that 

A(C) 



(3.15) 



and 



5(C) = -ln[l-MC)]- 



A(C) + 1 



InMC). 



(3.16) 



(3.17) 



1-MC) 

Fig.|| and Fig.H show A(C) and S{Q, respectively. 

Finally, we must integrate over p (or (). Note that, for a surface area A 
in configuration space, the density of modes in momentum space is A/ (27r)^. 
(Since the shape of the "horizon" is rather than S"^ in this model, A and 
the total entropy are infinite. However, we can perfectly define the entropy 
per unit area, and consequently we can pretend as if A is finite in our formula 
for S.) Therefore, 



S 



A 



(2vr)^ 

2tt Jo 

ttA 



ir/a 



d'p s{C) 

pdp S{C) 



2a2 Jo 



CdC 5(C), 



(3.18) 



where we have used the assumption that ma ~ mlpi <^ 1. Note that, even 
though S{() — > oo for C — ^ 0, CS{() becomes zero at C = 0. Thus, one finds 



where 



C=^ /'cdC 5(0 - 0.057. 

/ JO 



(3.19) 
(3.20) 
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If we consider that the quantum fluctuation of the horizon is 

a ~ 2v^/p; ~ 0.48 X /pi, (3.21) 

then the entanglement entropy is consistent with the Bekenstein-Hawking 
entropy Eq.( |l.l| ). 

3.2 More Realistic Model 

Next, we consider a free scalar field in the flat spacetime but adopt the 
Rindler coordinates {T,^,y,z) = {t,^,x), which are defined by 

t = ^ sinh ar, 

X = ^coshar, (3.22) 



where a is a constant [^]. The Rindler coordinates cover only a quarter of 
the Minkowski spacetime, x > \t\, called the Rindler wedge. The boundary of 
this Rindler wedge ^ = is the horizon for a uniformly accelerated observer 
in the Rindler wedge. In the Rindler coordinates, the flat metric becomes 

ds^ = -dt^ + dx^ + dy^ + dz^ = -i'^a^dr'^ + c/^^ + dy"^ + dz^ . (3.23) 

On the other hand, the most general, static and spherically symmetric 
black hole in four dimensions is 

ds"^ = -fir)dt^ + -r^c/r^ + R^ir) (dO^ + sin^ 9d(p'^), (3.24) 

where the horizon is at r = Vh which satisfies f{rh) = 0. We thus make the 
coordinate transformation [t, r, 6, (p) —>■ {t, rj, 6, 0), which is defined by 

^^1/7, dr, = Y^^dr, (3.25) 

where 

n = ldrf\r=r, (3.26) 

is the surface gravity of the black hole. Note that the horizon is at 77 = 0. 
The metric ( 3.24| ) becomes 



ds^ = -rfn" df + ^^^drf + R^{r) idO^ + sin^ Odcj)''), (3.27) 



16 



especially near the horizon 77 — 0, 

ds"^ -rfn^ dt^ + drf + R^{r) {dO'^ + sin^ 



(3.28) 



Therefore, by the comparison of Eq. (|3.23| ) with Eq.( |3.28D , we can think of 
the Rindler wedge as the model for the black hole, even though the shape of 
the horizon is now rather than S'^. 

Since we think of the Rindler time r as the "time", the index A is now 
replaced by (^,a;). The horizon is at ^ = and the thin region near the 
horizon becomes A = {(^,a;)|0 < ^ < a}. (See Fig.|^.) The "ground state" is 
the Rindler vacuum (rather than the Minkowski vacuum) which corresponds 
to the Killing vacuum in the case of a black hole. This is because the "time" 
is the Rindler time r rather than the ordinary Minkowski time t. 

Note that this horizon ^ = is a null surface, similar to the model of 
Ref. pSf. Thus, this model is not influenced by the criticism pO[ which is 



related to the fact that the boundary of the previous works like Refs. [|^, [19 
was timelike rather than null. 

The action of the scalar field in the Rindler coordinates becomes 



S = j drd^d^x -^a 



(3.29) 



Then, by using the orthogonality relations PB| , 



vr^ Jo 



°° dx 



X 



Ki^{x)Ki^{x) 



vr^ Jo 



du {21/ sinh Tru) Ki,^{x)Kii,{x') 



2z/sinh7rz/ 

X 6{x — x'), 



(3.30) 
(3.31) 



which are used in the Rindler quantization 1^^! o'^s finds that Gab and Vab 
appearing in Eq.( |2.14| ) become 



V{^,x;^',x') 



6{^-^')6{x~x'), 



e« 6i^-05ix-x'), 



(3.32) 

(3.33) 



a 



dv 



(2^ 



(2z/sinh7rz/) 



X K,,{M,i) KUM.O e^^ (^-^'). (3.34) 



17 



Then, one obtains that 



1 r du r d'^k 



{2v sinh ttz/) 



7 7r2 7 (27r)2 

X uKUMkO KUMkO e'^-(^-*'), (3.35) 
(iz/ r d^k 



W-\^,x;C,x') = /^/-|^(2.sinh.. 



7r2 y (27r)2 

X KUMkO KUMkC) e*^ (^-^'). (3.36) 



Thus, A"j, in Eq. (|2.42| ) becomes 
A(e,a^;r,^') = 



d^dx" [W-\^, x; e, x")W{e, x"; x') 

-6{^-0S{x-x'), (3.37) 



where (^, a;), (^', x') G A. 

By making the ansatz as above, 

Fi^,x)=e^P-^f{0, 
the eigenvalue equation reduces to 

f°- di' f"' dC [ dp , , ^ f dp' , , , 

^ /s:,„(A/,e) ifi^lAfp^") lUAMri") KiAMX) /(?') 

Then, by using the approximation, 

r diG{i) aG{a/2), 
Jo 

we find that 



(3.38) 



(A + l)/(0 (3.39) 



(3.40) 



. f dv , , . f dv' , , , , 
A + 1 = 4 / — (2z/smh7rz/) / — (2z/' smh ttz/' 



du' 



X - [ir,,(Mpa/2)]' [ir,,,(Mpa/2)]' . (3.41) 
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This would diverge unless we introduce a momentum cut-off for v and v' 
integrals. As in Eg. ( P.13 ), we decide the momentum cut-off in relation to 
the width of the region and the approximation Eq.( |3.40| ). By using ( p. 311 ), it 
can be decided by 



°- dv 
^ Jo vr 
''•''^ dv 



2 (2z/sinh7rz/)iri^(MpO Ki^{Mpa/2) 



~ 2/ — (2z/sinh7ri/) [J^i^(M„a/2)]^ (3.42) 

JO TX 

Unfortunately, this integral can not be done analytically. However, we can 
perform the numerical integration. Note that the cut-off is not a constant 
but a function of C = Mpo/vr by the dimensional analysis. Thus, 



A(C) 



^ f [MO dv . , , ["•^iO dv' , , 
4< / —{2vsmhTTVj / -^[^^ smhTrz/) 
[Jo 71 Jo TT^ 



X - [KU</2)f [K,,'inC/2)f \ - 1, (3.43) 



where t'c(C) is defined by 

rMO dv 9 1 

/ — (2z/sinh7rz/) [KU</2)f = -. (3.44) 

Jo TT"^ 2 

Then, from Eg. ( p.52| ) and Eg. (|2.53|) , one obtains /i(C) and S{Q, as above. 
A(C) and S{() are shown in Fig.|^ and Fig.^, respectively. (Note that X{() 
and S{() seem to be not smooth at C ~ 0.02 or 0.16. However, this is because 
the cut-off Vc{C)y which is shown in Fig.^ varies so rapidly there. Thus, we 
need more accuracy at such points.) 

Finally, after integrating over p (or () by using the fact that the density 
of modes in momentum space is A/{27cy, one finds that 

5'~c4, (3.45) 

where ^ 

C=-f Ct^C '5(C) ~ 0.089. (3.46) 
2 Jo 

(Even though A and 5* in our formula, if literally taken, are infinite since the 
shape of the horizon is rather than 5*^ in this model, we can still define 
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the entropy per unit area precisely. We can thus pretend as if A and S are 
finite in our final formula.) If we consider that the quantum fiuctuation of 
the horizon is 

a r-^ 2VC Ipi 0. 60 X Ipi, (3.47) 

then the entanglement entropy is consistent with the Bekenstein-Hawking 
entropy Eq.( |l.l|) . 

4 Conclusion and Discussion 

In this paper, we have considered the entanglement entropy between the 
outside and the thin region (of the order of the Planck length) of the inside 
the horizon based on the discussion in Sec.|l|. By constructing two models, 
a simple one and a more realistic one, we have shown that its entanglement 
entropy becomes 

5~c4, (4.1) 

where a is the quantum fiuctuation of the horizon and C is a constant. If 
the quantum fiuctuation of the horizon is 

ar^2VClpi, (4.2) 

we can interpret the Bekenstein-Hawking entropy, Eq. (|1 . 1|) , in the context of 
the entanglement entropy. This is consistent with the assumption that the 
quantum fiuctuation of the horizon is of the order of the Planck length. 

Although some authors have considered the entanglement entropy as the 
correction to the Bekenstein-Hawking entropy generated by matter fields, we 
want to consider this entanglement entropy as the Bekenstein-Hawking en- 
tropy itself. This is because we have considered the entanglement entropy 
of the Rindler vacuum (rather than the Minkowski vacuum) in the second 
model, which does not contain the thermal radiation of the Rindler particles. 
In the case of a black hole, this corresponds to the Killing vacuum (rather 
than the Kruskal vacuum), which does not contain the Hawking radiation. 
Thus, this entropy is not associated with the existence of the thermal radia- 
tion of particles but rather with the existence of the black hole itself. 

One might think that this entanglement entropy would depend on the 
number of matter fields which are present in the real world. That is, if there 
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are matter fields independently, one might think that the entanglement 
entropy would be multiplied by and conclude that this entropy could not 
be considered as the entropy of the "black hole" , since it would depend on A^. 
However, the entanglement entropy of the horizon in fact does not depend on 
N. This is because the quantum fluctuation of the horizon a also depends on 
and, besides, it is proportional to \/N. This can be seen from the following 
argument. Let us consider a Schwarzschild black hole with its mass M, which 
fluctuates within 6M {6M/M -C 1). Then, the Schwarzschild radius of this 
black hole fluctuates within 26M in the coordinate length. The proper length 
of this fluctuation becomes 

r=2{M+5M) rr=2{M+SM) rl^ 

ds= / , ^2^2M5M. (4.3) 

r=2M Jr=2M Wl-2M/r 



NxC^^ — = C— , (4.4) 



Note that 5M is proportional to A^, since the rate of spontaneous quantum 
emission or absorption of particles is proportional to A^. Thus, the fluctuation 
of the horizon is proportional to -y/iV in the proper length. (This is similar to 
the "brick wall" of 't Hooft |3^.) Since the coefficient in front of dS,'^ is 1 in 
Eq.( |3.23| ), a is the proper length and thus is proportional to \/N . Therefore, 
if the species of matter fields becomes A^, the entanglement entropy of the 
horizon becomes 

which is independent of A^ [Q. We thus consider this entanglement entropy 
of the horizon as the entropy of the "black hole" itself rather than the "matter 
field". 

The result of our analysis suggests that we can consider that the infor- 
mation on the collapsed star is stored as the EPR correlation between the 
outside and neighborhood (of the order of the Planck length) of the hori- 
zon. Since the horizon remains stable to the Planck scale, we can encode 
the enormous information on the collapsed star. If we used an ordinary wall, 
we could not do so because it begins to fluctuate far below the Planck scale. 
The information available outside the horizon is the probability distribution 
of the effective states (the effective density matrix) when we ignore the field 
near the horizon. Note that this consideration does not contradict with the 
no-hair theorem. 

Moreover, this picture appears to be consistent with the Euclidean path- 
integral approach by Gibbons and Hawking ^, |TT]. The entropy in that 
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approach arises from the fixed point of the Euchdean time translation or 
non-trivial topology of (r, r) section. In our analysis, we find that the quan- 
tities which appear in the first law of the black hole thermodynamics can be 
understood in relation to the Euclidean time translation: That is, the energy 
is its charge, the temperature is its period and the entropy is concerned with 
its fixed point. 

Finally, to be more realistic, we have to consider the case where the 
shape of the horizon is S'^, like Eg. ( p. 281) . We then have to expand the field 



by the spherical harmonics Yim{0,(f)). However, we expect that this would 
not change the result drastically and would turn out to be consistent with 
the result of Ref. p5[| , 

S ~ 0.024 X — . (4.5) 
a'' 

This is because as long as the radius of the sphere is much larger than the 
Planck length Ipi (which is equivalent to the near- horizon limit), we can ap- 
proximate the horizon as a plane. Of course, by using the method developed 
in this paper which is based on the Bombelli-Koul-Lee-Sorkin type calcula- 



tion |jT8| rather than the Srednicki type calculation |T9|, we will be able to 
obtain the result in a much simpler and more direct way. This is left for a 
future work. 
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Figure 3: The numerical evaluation for X{() 
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Figure 4: The numerical evaluation for S{Q 
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Figure 5: The numerical evaluation for VciC) 
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